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INTRODUCTION 


The  economic  or  military  value  of  any  vehicle  depends  fundamentally 
on  its  ability  to  traverse  a specific  path  between  its  point  of  depar- 
ture and  its  destination.  Means  for  the  control  of  the  path  vary  widely 
and  depend  on  a variety  of  constraints.  A train,  for  example,  is  con- 
strained to  move  along  a track  but  is  not  steered.  Its  control  freedom 
is  merely  one  of  speed.  An  automobile  or  surface  ship,  on  the  other 
hand,  while  constrained  to  move  on  the  surface  of  the  land  or  the  sea, 
must  be  steered  as  well.  Underwater  vehicles,  however,  have  an 
additional  degree  of  freedom;  for  this  reason  their  problems  of  the 
control  are  of  unusual  complexity. 

An  underwater  vehicle  or  weapon  system  contains  spatial  sensors, 
and  guidance  and  control  devices  (possibly  all  contained  in  the  human 
pilot)  whose  purpose  it  is  to  develop  three-dimensional  flight  path 
commands  appropriate  to  steering  co  as  to  reach  a destination  or  target, 
and  then  execute  those  commands  by  maintaining  or  modifying  the  forces 
on  the  vehicle  so  as  to  maintain  or  modify  the  velocity  vector.  This 
allows  the  intended  purpose  or  mission  to  be  accomplished.  Qualities 
of  a vehicle  that  tend  to  make  it  resist  changes  in  the  direction  or 
magnitude  of  its  velocity  vector  are  referred  to  as  stability,  while 
the  ease  and  expedition  with  which  the  vector  may  be  altered  are 
referred  to  as  control.  1 

The  purpose  of  this  report  is  to  develop  the  equations  of  motions 
and  their  associated  transfer  functions  for  underwater  vehicles  in  a 
form  suitable  for  use  by  engineers  charged  with  solving  these  complex 
problems  of  stability  and  control. 

One  begins  by  equating  the  forces  and  moments  (hydrodynamic, 
gravity,  and  buoyancy)  acting  on  the  vehicle  to  its  reactions,  in 
accordance  with  Newton's  laws.  The  body  is  assumed  to  be  rigid  with 
the  axis  system  located  at  its  c.g.  The  external  forces  and  moments 
acting  on  the  vehicle  are  represented  by  a Taylor  series  expansion 


1 McRuer,  Duane;  Ashkenas,  Irving;  and  Graham,  Dunston,  Aircraft 
Dynamics  and  Automatic  Control , pp.  203-220,  Princeton  University 
Press,  1973. 
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about  an  equilibrium  condition.  The  theory  of  small  perturbations  is 
introduced  into  the  derivation  to  arrive  at  linear  decoupled  equations. 
These  equations  can  be  solved  in  either  the  time  domain  2 or  in  the 
frequency  domain.  The  technique  of  frequency  domain  solution  using 
Laplace  transforms  is  chosen  here.  It  is  then  possible  to  use  conve- 
nient transfer  function  models  for  the  dynamics  of  the  vehicle  and  all 
the  analytical  techniques  for  the  study  of  linear  feedback  systems  can 
be  brought  to  bear  on  the  problem.  Appropriate  nondimensional  factors 
are  introduced  to  conform  to  the  standard  nomenclature  followed  in  Navy 
publications.  3 The  complete  six  degree  of  freedom,  small  perturba- 
tion, linearized,  decoupled  equations  of  motion  are  then  presented  for 
a self-propelled  vehicle.  An  example  is  given  using  a typical  under- 
water vehicle  to  demonstrate  the  application  of  the  equations. 


INERTIAL  FORCES  AND  MOMENTS 


The  equations  of  motion  employed  are  a linearization  of  the 
classical  Euler  equations  for  a rigid  body  with  respect  to  a set  of 
axes  fixed  at  the  body  c.g.  The  positive  directions  of  the  axes, 
angles,  linear  velocity  components,  angular  velocity  components, 
forces,  and  moments  are  shown  in  Figure  1.  Since  the  derivation  of 
these  equations  can  be  found  in  any  standard  text  on  dynamics,  the 
details  will  not  be  presented  here.  The  reader  should  refer  to  McRuer, 
et  al.  for  a complete  derivation  of  the  inertial  forces  and  moments. 
Referred  to  body-fixed  axes  one  may  write 

LX  = m(U  + QW  - RV) 

EY  = m(V  + RU  - PW)  (1) 

'Ll  = m(W  + PV  - QU) 


2 Naval  Ship  Research  Development  Center,  Report  P-433-M-01,  Users 
Guide,  NSRDC  Digital  Program  for  Simulating  Submarine  Motion  ZZMN- 
Revision  1.0,  by  Ronald  W.  Richards,  June  1971. 

3 The  Society  of  Naval  Architects  and  Marine  Engineers,  TMB  No.  1-5, 
Nomenclature  for  Treating  the  Motions  of  a Submerged  Body  Through  a 
Fluid,  April  1952. 

1 ibid. 
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The  following  assumptions  were  made  in  the  derivation  of  these 
equations : 


1.  The  vehicle  is  assumed  to  be  a rigid  body. 


2.  The  mass  and  mass  distribution  of  the  vehicle  are  assumed  to 
be  constant. 


3.  The  XZ  plane  is  assumed  to  be  a plane  of  symmetry. 
These  equations  are  programmed  in  Reference  2. 


TRANSFORMATION  MATRIX 


For  application  of  the  equations  of  motion  to  underwater  vehicles, 
the  effects  of  earth  rotation  can  be  neglected  and  the  earth  can  be 
assumed  to  be  flat.  These  assumptions  allow  one  to  consider  the 
initial  local  level  coordinate  frame  (x  pointed  north,  y pointed  east, 


o o 

and  z pointed  down)  to  be  an  inertial  coordinate  frame.  The  gravity 

forces  are  known  in  the  initial  local  level  coordinate  frame  and  the 


hydrodynamic  forces  are  measured  in  a body  coordinate  frame.  To  take 
advantage  of  available  hydrodynamic  test  data,  the  equations  of  motion 
will  be  referred  to  the  body  coordinate  frame.  Therefore,  the  forces 
measured  in  the  initial  local  level  coordinate  frame,  i.e.  gravity, 
must  be  transformed  to  the  body  frame.  This  is  accomplished  by  a 
series  of  rotations  through  Eulerian  angles  as  shown  in  any  standard 
text  on  dynamics.  The  transformation  matrix  which  performs  this  is 
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LINEARIZATION  OF  THE  EQUATIONS  OF  MOTION 


Equations  (1)  are  to  be  linearized  before  they  are  expanded  to 
include  the  hydrodynamic  and  gravity  forces.  These  equations  contain 
products  of  the  dependent  variables,  therefore  they  are  in  general 
nonlinear.  To  reduce  them  to  a tractable  form,  the  total  motion  can 
be  considered  as  composed  of  two  parts:  an  average  or  mean  motion 

that  is  representative  of  the  operating  point  or  equilibriums  con- 
ditions, and  a dynamic  motion  that  accounts  for  small  perturbations 
about  the  mean  motion.  Accordingly,  each  of  the  total  instantaneous 
velocity  components  of  the  vehicle  can  be  written  as  the  sum  of  a 
velocity  component  during  the  equilibrium  condition  and  a change  in 
velocity  caused  by  the  disturbance: 

U = U + u 
o 

V = V + v 
o 

W = Wo  + w (3) 

P = P + p 
o 

Q = Qo  + q 

R = R + r 
o 

The  zero  subscripts  in  Equations  (3)  indicate  the  steady  flight 
velocities,  and  the  lower  case  letters  represent  the  changes  in  the 
velocities  (disturbance  velocities). 

During  the  equilibrium  condition,  the  vehicle  is  assumed  to  be 
flying  with  wings  level  and  with  all  the  components  of  velocity  zero 
except  Uq.  Thus,  Equation  (3)  may  be  rewritten  as 

U = U + u 

o 


V = v 
W = w 
P = P 
Q = q 
R = r 


(4) 
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Substituting  Equations  (4)  in  Equations  (1)  yields 
IX  » m(u  + qw  - rv) 


IY  = m(v  + rU  + ru  - pw) 

IZ  = m(w  + pv  - qU^  - qu) 

ZK  = pi  - rl  + qr(I  - I ) - pql 

x xz  ^ z y xz 

EM  = ql  + pr (I  - I ) - r2I  + p2I 

y x z xz  xz 


ZN  = rl  - pi  + pq(I  - I ) + qrl  . 
z xz  r y x xz 


The  disturbances  from  the  equilibrium  flight  condition  will  now  be 
assumed  to  be  small  enough  so  that  the  products  and  squares  of  the 
changes  in  velocities  are  negligible  in  comparison  with  the  changes 
themselves.  Also,  the  disturbance  angles  are  assumed  to  be  small 
enough  so  that  the  sines  of  these  angles  may  be  set  equal  to  the  angles 
and  the  cosines  set  equal  to  1.  Products  of  these  angles  are  also 
approximately  zero  and  can  be  neglected.  Thus,  terms  similar  to  rv  and 
uq  may  be  set  equal  to  zero,  and  Equations  (5)  reduced  to 

ZX  = mu 

ZY  = m(v  + rU  ) 
o 

ZZ  = m(w  — qU  ) 

O / r \ 


£K  = pi  - rl 

r x xz 

ZM  = ql 


ZN  = rl  - pi  . 
z xz 


The  relationship  between  the  angular  velocities  and  the  rate  of 
change  of  the  angles  ($,0,\j/)  are  given  by 

P = $ - i sin  8 

Q = 0 cos  <f>  + ip  cos  0 sin  ij> 


R = cos  0 cos  <j)  - 0 sin 
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body  axis  and  <f> , 0,  and  ij;  are  the  angular  velocities  measured  about  the 
equilibrium  axis  system.  P,  Q,  and  R are  the  angular  velocities 
1 measured  by  a rate  gyro  located  in  the  vehicle.  Applying  Equations  (4) 

and  the  small  perturbation  assumption,  these  equations  reduce  to 

p = $ 

q = 9 (7) 

r = 

The  small  perturbations  assumption  not  only  limits  the  applica- 
bility of  Equations  (6)  to  small  changes  in  angles  and  velocities  about 
an  equilibrium  condition,  but  reduces  Equations  (6)  to  linear  equations 
and  yields  a simplification  of  the  mathematical  methods  necessary  for 
the  analysis  of  the  complicated  vehicle  motions.  In  a rigorous  mathe- 
matical sense.  Equations  (6)  are  applicable  only  to  infinitesimal  dis- 
turbances; however,  experience  has  shown  that  quick  and  accurate  results 
can  be  obtained  by  applying  these  equations  to  disturbances  of  more 
finite  magnitude. 

EXPANSION  OF  THE  HYDRODYNAMIC  FORCES  AND  MOMENTS 


The  hydrodynamic  forces  are  exerted  on  the  vehicle  by  the  surround- 
ing fluid.  They  are  present  whenever  there  are  any  reactive  forces 
between  the  fluid  mass  and  the  vehicle.  In  steady  flight  they  result 
from  relative  motion  between  the  vehicle  and  the  fluid  mass  or  from 
accelerated  flows  produced  by  the  deflection  of  a control  surface. 
Although  the  specific  forces  depend  on  their  peculiar  origins,  the  form 
of  the  expressions  that  describe  perturbations  in  these  forces  is  not 
particularly  dependent  upon  origin.  For  example,  on  vehicles  designed 
to  be  propelled  at  very  low  speeds  by  divers,  the  dominant  forces  and 
moments  are  produced  by  the  accelerated  flows  (added  mass  and  moments 
of  inertia)  and  the  gravity-buoyancy  terms  (metacentric  moment) . Where- 
as in  the  case  of  a high  speed  submarine,  the  dominant  contributions 
come  from  the  steady  hydrodynamic  terms.  The  end  results  of  the  treat- 
ment here  are  pertinent  to  all  kinds  of  underwater  vehicles. 

These  forces  and  moments  are  known  to  be  functions  of  the  relative 
velocity,  acceleration,  and  position  as  well  as  control  deflections  1 . 
In  functional  form  they  are  written  as 


1 ibid. 
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f(u,  v,  w,  u,  v,  w,  q,  q,  0,  p,  p. 


r. 
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If  the  hydrodynamic  forces  are  considered  to  be  continuous  functions  of 
all  these  variables,  each  of  the  forces  (X^,  Yj  , and  Z^)  and  the  moments 
(K^,  M^,  and  N^)  can  be  expressed  in  terms  of  the  variables  by  a Taylor 
series  expansion  about  the  equilibrium  condition.  Because  of  the  small 
perturbation  assumptions,  second  order  and  higher  terms  of  the  Taylor 
series  are  neglected.  Additionally,  because  the  XZ  plane  is  a plane  of 
symmetry,  X^,  Z^,  and  M^  are  functions  of  only  u,  w,  q,  their  deriva- 
tives and  0,  whereas  , and  are  functions  of  only  v,  p,  r,  their 

derivatives  and  if.  Thus,  expanding  the  forces  and  moments  about  an 
equilibrium  point  yields 
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Y , Zq,  Kq,  Mq,  and  NQ  are  the  forces  and  moments  acting  along  and 
about  the  X,  Y,  and  Z axes,  respectively,  while  the  vehicle  is  in  the 
equilibrium  flight  condition.  The  terms  similar  to  Xuu  express  the 
change  in  the  given  force  and  moment  caused  by  the  disturbance  quantity. 
The  term  Xu  is  known  as  either  a stability  derivative  or  a hydrodynamic 
coefficient  and  is  defined  as  the  change  in  the  X-force  with  respect  to 
the  u-velocity,  evaluated  at  the  equilibrium  condition.  That  is, 

X = (f ) . 
u 3u  o 


EXPANSION  OF  THE  GRAVITY  AND  BUOYANCY  FORCES  AND  MOMENTS 


The  gravity  force  can  be  considered  to  act  at  the  center  of  gravity, 
whereas  the  buoyant  force  acts  at  the  center  of  buoyancy.  The  buoyant 
force  also  creates  a moment  since  it  is  displaced  from  the  c.g.  by  the 
distance  X„,  Yg,  Zg.  These  distances  are  shown  in  Figure  2 for  a stable 
and  unstable  configuration.  To  find  the  expressions  for  the  components 
of  gravity  and  buoyancy  to  be  used  in  the  equations  of  motion,  the 

(Text  Continued  on  Page  10) 
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general  transformation  matrix  is  again  used  to  transform  these  forces 
from  the  initial  equilibrium  flight  conditions  to  the  disturbed  flight 
condition. 


The  components  of  force  due  to  gravity  and  buoyancy  acting  on  a 
vehicle  which  is  in  equilibrium  flight  and  which  has  initial  angles  6 
and  d>0  with  respect  to  the  gravity  vector  can  be  found  by  direct  reso- 
lution of  the  force  along  the  equilibrium  flight  axes  Xq,  Yq,  Zq: 

X = - (W  - B)  sin  9 
oGB  o 


Y = (W  - B)  cos  0 sin  $ (9) 

oGB  o o 

Z = (W  - B)  cos  0 cos  <J>  . 

oGB  o o 

The  components  of  force  due  to  buoyancy  alone  are 
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(10) 
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- B cos  0 
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The  components  of  gravity  and  buoyancy  which  act  along  the  disturbed 
axes  x,  y,  and  z can  be  found  by  substituting  Equations  (9  and  10)  into 
the  transformation  matrix  given  in  Equations  (2).  Thus, 
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or 


XGB  = XoGB  COS  * cos  6 + YoGB  Sln  * C0S  6 " ZoGB  Sln  0 


V = X (cos  <P  sin  0 sin  <p  - sin  ip  cos  $) 

+ yoGB(cos  ^ cos  Q + sin  <P  sin  6 sin  <J>)  + Z cos  6 sin 


(11) 


oGB 


(cos  ^ sin  0 cos  <(>  + sin  ip  sin  (f>) 

+ YQ(,^(sin  0 cos  $ sin  - sin  <f>  cos  4/)  + Z ^ cos  0 cos  <f> 
Applying  the  small  perturbation  assumption,  these  equations  reduce  to 

XGB  = XoGB  + ^YoGB  ~ 0ZoGB 


YGB  “ ^XoGB  + YoGB  + <f>ZoGB 


(12) 


JGB 


0XoGB  ~ <(>YoGB  + ZoGB  * 


In  a similar  fashion,  the  moment  equations  are  found  from 

XB  " XoB  + *YoB  " 9ZoB 

YB  " - *XoB  + YoB  + *ZoB 

Z_  = 0X  t - <j>Y  D + Z _ 

B oB  oB  oB 

Md  = zc  - ZD  xD  = z„  (X  D + iJjY  d - 0Z  „)  - xn(0X  _ - <f>Y  _ + Z D) 
B B B B B B oB  oB  oB  B oB  oB  oB 


(13) 


NB  = "XB  yB  + YB  XB  = -yB(XoB  + *YoB  " 0ZoB)  + V'  *XoB  + YoB  + ♦W 
= ZR  yn  - Y„  = yR(0X„n  - *Y  n + Z - zn(-  + Y.n  + *Z  • 


Since  it  is  assumed  that  the  wings  are  level  in  equilibrium  flight, 

<P  = 0 and  thus  YqB  will  be  zero.  However,  0 is  arbitrary  since  a 
vehicle  can  climb°or  dive  in  steady  flight,  furthermore,  since  the 
XZ  plane  is  assumed  to  be  a plane  of  symmetry,  Y = 0.  Thus,  Equations 
(12)  and  (13)  reduce  to 
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*CB  ' XoCB  ' 8<”  - “>  “o 


Y„D  » + iHW  - B)  sin  0 + 4> (W  - B)  cos  0 

Gd  o o 


ZCB  « sl”  °o  + ZoGB 

(14) 

M=Xz-xZ+  0z  B cos  0 - 0x  B sin  0 

BoBBBoB  B o B o 

Nu  = - x_B<J>sin  0 - xuB<)>cos  0 

D D O D O 

Kd  = z„i|/Bsin  0 + z <}>Bcos  0 

B B o B o 

It  is  conventional  in  Navy  literature  3 to  treat  the  gravity  and 
buoyancy  contributions  as  hydrodynamic  coefficients.  Thus, 

XGB  “ XoGB  + V 

YGB  = YoGB  + V + V 


ZGB  ZoGB  + Z09 


MB  =MoB+M00 
NB  = NoB  + V + V 


KB  = KoB  + V + V 


Comparison  of  Equations  (13)  and  (14)  shows  that 

XQ  = - (W  - B)  cos  0 


Y = (W  - B)  cos  0 
<p  o 

Y,  = (W  - B)  sin  0 
Za  = - (W  - B)  sin  0 

D O 

Mq  = Bzd  cos  0 - Bxn  sin  0 

y d o d c 


3 ibid. 
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N = - x„B  cos  0 
<J>  B o 

N = - x B sin  0 

* ‘ ° (16) 

K = z B cos  0 Con. 

B o 

K = z_B  sin  0 
^ B o 

Therefore,  in  the  following  section  when  the  contributions  are 
grouped  together  to  form  the  complete  equations  of  motion,  the  gravity 
and  buoyancy  disturbance  terms  XQ0,  Y d>,  etc,  will  be  included  in  the 
hydrodynamic  disturbance  terms.  ^ 


COMPLETE  EQUATIONS  OF  MOTION 


The  individual  contributions  to  the  equations  of  motion  have  been 
examined  in  detail.  The  complete  equations  of  motion  of  the  vehicle 
can  now  be  written.  This  is  accomplished  by  summing  the  contributions 
of  the  external  forces  (Equations  8 and  15)  and  equating  this  sum  to 
the  right-hand  side  of  Equation  (6) . 

The  equations  for  the  equilibrium  flight  condition  can  be  found 
by  substituting  the  equilibrium  values  of  the  hydrodynamic  and  gravity 
forces  and  moments  into  Equations  (6)  and  setting  the  disturbance  terms 
equal  to  zero: 


X + 
o 

Y + 
o 

Z + 

o 

K.  + 
o 

M + 
o 

N + 

o 


XoGB  ° 


YoGB  " ° 


ZoGB  " ° 


K n = 0 
oB 


M _ = 0 
oB 


N D = 0 
oB 


(17) 


The  equations  of  motion  for  the  disturbed  vehicle  are  found  by 
substituting  tin  disturbed  values  of  the  forces  and  moments  Equations 
(8)  and  (15),  into  Equations  (6): 
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+ (X1i) 


Hydrodynamic 


oGB 


+ (X.i) 


Gravity-Buoyancy 


The  quantities  in  the  boxes  sura  to  zero  because  of  the  equilibrium 
flight  condition  of  Equations  (17).  Hence, 


I 


mu 

= (Xti) 

+ (X  i) 
H 

GB 

Z,  K, 

M,  and 

N equations 

• 

ra(w  - 

U 0)  = 
o 

(Z  i)  + 
H 

(Z.i) 

I 0 = 

y 

(Mi)  + 
H 

(M.i) 

ra(v  + 

v>  ■ 

(Y  i)  + 
H 

(Y.i) 

I 4 - 

X 

I H . 

XZ 

(K  i)  + 
H 

(KLi) 

I T - 
z 

14  = 

XZT 

(N  i)  + 
H 

(N.i) 

GB 


GB 


GB 


GB 


GB 

where  the  terras  (X^i)  represent  the  hydrodynamic  coefficients  of 
equations  (8)  and  (X^i)  represent  the  gravity-buoyancy  coefficients 
of  equations  (15).  With  the  force  and  moment  terms  expanded  the 
equations  become 

mil  = X-u  + X u + X»w  + X w + X-q  + X q + Xn0  + Xx  6 
u u w w q q 0 os 

s 

m(w  - II  0)  = Z*u  + Z u + Z*w  + Z w + Z*q  + Z q + Z.0  + Zx  5 (18) 

o u u w w qn  qn  0 6^  s 

10=  M* u + M u + M*w  + M w + M*q  + M q + M O + Mr  6 
y u u w w q q 0 °s  S 

and 
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L 


ra(v  + U 4))  = Y*  v + Y v + Y* p + Y p + Y^  + Y-r  + Y r + Yr  6 
o v v p pr  4>  r r 6D  R 

K 

1$  ~ I * = K-v  + K v + K;p  + K p + K 4»  + K;r  + K r + K 6 (19) 

x xz  v v p p <p  r r uD  k 

k 

If  - I <f>  = N-v  + N v + N-p  + N p + N <J>  + N-r  + N r + N„  5„ 

2 xz¥  V V p*^  pK  <f>¥  r r oR  R 

Equations  (18  and  19)  are  the  linearized  small  perturbation  equations 
of  motion. 

Examination  of  these  equations  shows  that  Equations  (18)  are 
functions  of  the  variables  u,  9,  and  w,  whereas  Equations  (19)  are 
functions  of  the  variables  v,  r,  and  p,  thus  as  a result  of  the 
assumptions  made  in  the  previous  analysis,  the  equations  of  motion  can 
be  treated  as  two  independent  sets  of  three  equations.  Equations  (18) 
are  referred  to  as  the  longitudinal  or  symmetrical  equations  because, 
when  these  motions  occur,  the  plane  of  symmetry  of  the  vehicle  remains 
in  the  plane  it  occupied  in  the  equilibrium  flight  conditions.  Equa- 
tions (19)  are  referred  to  as  the  lateral  or  asymmetrical  equations. 
Since  the  longitudinal  motions  are  independent  of  the  lateral  motions, 
they  are  treated  separately  in  the  remainder  of  this  report. 


DERIVATION  OF  THE  TRANSFER  FUNCTIONS 


In  this  section,  Equations  (18)  and  (19)  are  converted  by  the  use 
of  determinants  into  transfer  functions.  The  Laplace  transform  method 
of  solution  is  used.  The  transformed  equations  are  nondimensionalized 
to  conform  to  the  standard  SNAME  nondimensionalizing  factors  given  in 
Reference  3.  Table  1 shows  a representative  list  of  these  factors.  It 
should  be  noted  that  the  resulting  transfer  functions  are  nondimensional . 
Before  evaluating  the  characteristic  frequency  or  the  time  response  of 
the  vehicle,  the  transfer  function  must  be  first  converted  to  the 
dimensional  form.  The  longitudinal  transfer  functions  are  derived 
first,  followed  by  the  lateral  transfer  functions. 

Longitudinal  Transfer  Functions 


Substituting  Equations  (7)  into  Equations  (18)  yields 


mu  - 

X«u  - 

X u - 

X*w  - 

X w - 

x*e  - 

X 0 - 

xQ0  = : 

u 

u 

w 

w 

q 

q 

0 

• 

• 

• 

• 

•• 

• 

• 

raw  - 

mU  e 

- Z-u 

- Z u 

- Z*w 

- Z w 

- Z*0 

- Z 0 - 

o 

u 

u 

w 

w 

q 

q 

I 0 - M*u 

-Mu 

- M*w 

- M w 

- M*  0 

-Me 

- M.0  = 

• 

y 

u 

u 

w 

w 

q 

q 

0 

I x u 

6 s 
s 


(20) 
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TABLE  1 

SNAKE  NOMENCLATURE  FOR  NONDIMENS IONAL  COEFFICIENTS 


m ’ 

m/^jJ l 3 

X’ 

w 

X /bpi2 u 
w o 

M’- 

w 

= 

M*  /^p 
w 

s’ 

= 

CD 

c: 

o 

X’-  = 

q 

x-Zhpi1* 

q 

M’ 

w 

= 

M /^f.3U 
w o 

u’ 

• 

u/U 

o 

X’ 

q 

\'*iX 

M'- 

q 

= 

w* 

w/U 

o 

x0  " 

x />spfc2u2 

b o 

M’ 

q 

= 

M /%pXAU 

q o 

X’* 

u 

- 

x^/'-spi3 

X*«  = 
s 

X.  /hpl2U2 
o O 

S 

Mo 

= 

H /hps>  3 V2 
b o 

X* 

u 

- 

x /hpa2v 

U O 

M * • = 

u 

M,5 

s 

= 

M.  /% pi3u 
c < 

s 

X’> 

w 

= 

X*/*jpJ.3 

W 

M* 

u 

M /‘2pf.3U 
u o 

I' 

y 

= 

i fhpi5 

The  right-hand  side  of  Equations  (20)  are  the  control  forces  and 
represent  the  means  by  which  either  a human  operator,  an  autopilot,  or 
a disturbance  input  can  control  the  motion  of  the  vehicle.  The  control 
surface  and  disturbance  inputs  are  the  forcing  functions  which  deter- 
mine the  resultant  motion  of  the  vehicle.  Since  the  vehicle  equations 
of  motion  are  linear  equations,  the  principle  of  superposition  may  be 
used  to  obtain  a solution.  For  example,  the  response  to  simultaneous 
application  of  sternplane  and  bowplane  deflections  can  be  determined 
by  calculating  the  response  to  each  of  these  deflections  separately  and 
then  adding  together  the  results  to  arrive  at  a complete  solution. 

In  this  report,  only  the  longitudinal  response  to  sternplane 
deflections  and  the  lateral  response  to  rudder  deflections  are  given 
detailed  analysis.  It  should  be  emphasized,  however,  that  the  mathe- 
matical techniques  of  solution  for  other  control  inputs  are  identical. 

Applying  the  Laplace  transform  to  Equations  (20)  yields 

[(m  - X*)s  - X ]u  + [-  X*s  - X ]w  + [-  X*s2  - X s - XJO  = X*  6 
1 u u w wJ  1 q q 0J  6 s 

s 

[-  Z*  s - Z ]u  t-  f(tn  - z*)s  - Z ]w  + t-  Z’S2  - (Z  + raU  )s  - ZJ9  = Zx  6 
1 u u w w q q o 0 5s 

s 


[-  M*  s - M ]u  + [-  M*  s - M ]w  + [(I 


M-)s  -Ms-  MJ0  = Mr  5 (21) 

q q 0 5 s 
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Equations  (21)  are  nondimensionalized  by  dividing  the  force  equa- 
tions by  -spJl^U2  and  the  moment  equations  by  ‘-spf.3U2.  Dividing  the  first 
bracket  in  the  X-force  equations  by  *sp22U2  yields 


[ (m  - X* ) s - X ] u X* 

u u , , m u . s 

l \ 'i  “* 


2 2 
U 

o 


u 

■*  T 


. 3 , 0 3 U /*  i n2„  JU 

'■ip  i hpi  o hpi-  u ■ 


= [(m*  - X'Os'  - X’Ju' 


Likewise , for  the  second  and  third  brackets  and  the  right-hand  side  we 

obtain: 


X* s - X ]w 
w w 


'spsiV 

o 


[-  x-s  - x s - xje 
q q e 

2 2 

%pl  U 

o 


X- 

w 

S 

X 

w 

'-ipi3 

u H 

o 

’-spi2U 

- X' • s ' 
w 

1 - X ' 

w 

]w' 

X* 

q 

2 

s 

X 

q 

hpz1* 

u2/*2 

o 

»ip23l 

- X'  • s 

q 

,2  _ x, 

s’  - : 

q 

e 


2 2 

ru 

o 


]e 


Finally  for  the  right-hand  side  of  the  equation 


Xx  6 
6 s 
s 


= X!  6.  . 


. „ 2tr2  “5  s 
%p2  U s 

o 


Therefore,  the  nondimensional  X-force  equation  becomes 

[(»’  - X'.)s'  - X;ju’  + [-  X’-s’  - X’]w’  + [-X'.s-2  - X’s’  - X’]6  =x’  «£ 


Similarly,  the  remaining  longitudinal  equations  can  be  placed  in  non- 
dimensional  form  to  yield: 

[-  z'js'  - z;iu'  + [(».'  - z';)s'  - z;N' 

+ [-  z'-s'2  - <z;  +«')='  - Z'e)9  - z;  6S 
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[-  M'  • s f 
u 


M? ]u'  + [-  M* • s 1 - M' ] w * 
u w wJ 


+ [(I'  - M' • ) s ' 2 - M's’  - M* ] 0 = Ml  6 

y q q e 6S  s 

These  equations  are  shown  in  Table  2 for  reference  purposes.  The  equa- 
tions predict  the  longitudinal  small  perturbation  motions  of  self- 
propelled  vehicles. 


TABLE  2 

NONDIMEN SI ONAL  LONGITUDINAL  EQUATIONS  OF  MOTION 

[(m’  - X ' • ) s ’ - X’]u'  - [X'-s'  + X'Jw’  - [X’-s’2  + X's’  + X']0  = X'  6 
1 ' u u w w q q 0 °s  S 

l-Z’.s’  - z;]u*  + [(m*  - Z’.)s'  - z;jw* 

- [Z-.s-2  + (Z-  + m’)s’  + Z']0  = Z'^ 

[-  M'-s'  - M']u’  - [M'-s'  + M;]w'  + [(!'  - M'-)s'2  - M^s ' - M’]0  = M> 


The  transfer  functions,  for  a given  input,  are  obtained  by  solving 
the  transformed  simultaneous  equations  of  motion  for  the  output  variable 
of  interest  with  all  other  inputs  considered  to  be  zero.  For  example, 
using  determinants,  we  can  directly  write  the  pitch  attitude  to  control 
input  transfer  functions  as 


e(s)  _ 

5g(s)  - 


(m’  - X 

* • ) s * 
u 

- X1 
u 

X’ 

' • s f 

w 

- X* 
w 

X6 

s 

- Z ' • s ' 
u 

- Z* 
u 

(m’ 

- Z’ 

'•)s' 

w 

- z’ 

w 

7 1 

s 

- M'  • s ' 
u 

- M' 
u 

- 

M?  s * - 
w 

- M’ 
w 

s 

if  - X1- 
a 

)s'  - 

X'  -X ' • s ’ - X' 
u w w 

- X'- 

q 

s'2  - X's 

q 

' ~ xe) 

Zf  *sf  - 
u 

Z' 

u 

(m1  - Z'-Js' 
N w 

- 

Z* 

w 

-Z*  s ' 

q 

- (Z'  + 

q 

m'  )s' 

+ ~'i 

M*  •s*  - 
u 

M' 

u 

- Mf • s * - M' 
w w 

(i'  - 

y 

M' • ) s ' 2 

q 

- M's' 

q 
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By  expanding  the  determinants,  the  transfer  function  can  be 
expressed  as  the  ratio  of  a numerator  polynomial  in  s over  a denomina- 
tor polynomial.  The  denominator  polynomial  A(s),  is  common  to  all  the 
longitudinal  transfer  functions  and  its  factors  determine  the  frequency 
and  damping,  or  time  constants,  of  the  individual  modes  of  motion.  The 
numerator  polynomials  depend  on  the  output  quantity  of  interest.  The 
general  polynomial  forms  of  the  primary  longitudinal  transfer  functions 
are  given  in  Table  3.  The  expressions  for  the  various  A,  B,  C,  etc. 
coefficients  of  the  equations  in  Table  3 is  given  in  Appendix  A in 
terms  of  hydrodynamic  coefficients. 


TABLE  3 

LONGITUDINAL  TRANSFER  FUNCTIONS 

.0 


N , 


0/6 


A0s’2  + Bes'  + C0 


^Long  As'^  + Bs'3  + Cs'2  + Ds'  + E 


w'/6  = — 

s A 


A s 
w 


3 2 

' + B s'  + 


C s'  + D 
w w 


Long  As'^  + Bs'3  + Cs'2  + Ds'  + E 


A s'3  + B s 


<2  + 


u'/6 


C s'  + D 

U *1 


S ^Long  As ' ^ + Bs'3  + Cs'2  + Ds'  + E 


Lateral  Transfer  Functions 


The  lateral  transfer  functions  are  derived  in  a manner  similar  to 
the  longitudinal  transfer  functions.  Substituting  Equations  (7)  into 
Equations  (19)  yields 


mv  + 

mU  - 

Y-v  - 

Y v - 

Y- <f>  - 

Y <f>  - 

' YJ  * 

Y-Y  - 

Y ip  = 

Y.  6 

0 

V 

V 

P 

P 

<t> 

r 

r 

6r 

R 

I <(>  - 

- I T 

- K*v 

- K v 

- K-<f> 

- K i 

- M 

- k«v 

- K * 

= Kx 

6p  (25> 

X 

XZ 

V 

V 

P 

p 

r 

r 

6r 

R 

I ¥ ■ 

- I (f> 

- N.v 

- N v 

- N*  <p 

- N 4> 

- N 4> 

- N-Y 

- N } 

= N. 

6d 

z 

XZ 

V 

V 

P 

P 

r 

r 

6r 

R 
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The  nond imensional , lateral  equations  of  motion  are  obtained  by 
applying  the  Laplace  transform  to  Equations  (25)  and  dividing  the  force 
equations  by  SsP K. 3U2  and  the  moment  equations  by  *5p£2U2.  Table  4 shows 
the  lateral  equations  in  terms  of  0,  4> , and  ip  where  0 i v/U  . 


The  lateral  transfer  functions  are  found  from  the  equations  of 
motion  by  expanding  the  numerator  and  denominator  determinants  as  pre- 
viously noted.  The  general  polynomial  forms  of  the  primary  lateral 
transfer  functions  are  given  in  Table  5.  Expressions  for  the  coeffi- 
cients A,  B,  C,  etc.  are  given  in  Appendix  B in  terms  of  hydrodynamic 
coefficients. 


TABLE  4 


NOND IMENSIONAL  LATERAL  EQUATIONS  OF  MOTION 

[(m'  - Y ' • ) s ’ - Y ' ] 0 + l-  Y'-s'2  - Y ' s ' - Y ' ] tp 

1 v v p p <p 

+ [-  Y'-s'2  + ( n'  - Y')s']iJ/  = Y'  6 
r r <5  R 

K 

[-  K'-s'  - K']0  + [(I'  - K'-)s'2  - K's’  - K']$ 

1 v v x p p <j> 


+ K-  Kz  ~ K’r)s’  -k;s']*-k:6r 

K 


[-  N'-s’  - N ' ] 0 + [(-  I'  - N 1 • ) s ' 2 - N ' s ' - N'] 
1 v v xz  p p <p 


+ [(i;  - N ' ^) s ' 2 - N^s']*  = ^ 6r 

R 
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TABLE  5 


LATERAL  TRANSFER  FUNCTIONS 


6/<5d  = 


s'(Ags'3  + BgS'2  + CgS ' + Dg) 


^ ^L at  s' (As' 4 + Bs ' 3 + Cs ' 2 + Ds ' + 


E) 


*/6„  = 


s’  (A , s ' 2 + B,s'  + C ) 
_ <t>  <t> 

, 4 _ .3  . „ .2 


^Lat  s' (As' 4 + Bs'3  + Cs'^  + Ds ' + E) 


,P/^D  = 


_R_ _jj_ 


A, s'3  + B s'2  + C , s ' + D 


<P  ip 


R A 


Lat  s' (As' 4 + Bs'3  + Cs ' 2 + Ds'  + E) 


EXAMPLE 

To  illustrate  the  use  of  the  equations  of  motion  and  transfer 
functions,  the  dynamic  characteristics  of  a typical  vehicle  are  evalu- 
ated. The  vehicle  length  is  49.33  feet  while  the  velocity  is  10.13 
feec  per  second.  The  hydrodynamic  coefficients  of  the  vehicle  are 
given  in  Table  6. 


The  transfer  functions  in  Tables  3 and  5 are  evaluated  by  substi- 
tuting the  numerical  values  of  the  hydrodynamic  coefficients  from  Table 
6 into  the  expressions  in  Appendix  A and  B.  The  numerical  values  for 
each  transfer  function  coefficient  are  given  in  Table  7.  These  are  the 
nondimensional  transfer  function  coefficients.  To  obtain  the  dimensional 
form  of  the  transfer  function,  substitute  s'  = s i/U  into  the  equations 
in  Tables  3 and  5.  For  the  velocity  transfer  functions  u'/<5g,  v'/<$^, 
and  w'/<5  , the  left  side  of  the  equations  in  Tables  3 and  5 must  also 
be  multiplied  by  U to  convert  to  the  dimensional  form.  For  illustra- 
tive purposes  the  nondimensional  transfer  functions  8/6g  and  v'/6  will 
be  converted  to  the  dimensional  form.  The  nondimensional  6/6 g transfer 
function, 


Q/J.  1. 10~4( . 3291  s'2  + .4746  s'  + .1361) 

s - 1.10  (.0902s1  + .3960s*  + .3288s'  + .1028s'  + .0115) 

is  converted  to  the  dimensional  form  given  below  by  the  substitution 
s'  = s A/U  . 


(Text  Continued  on  Page  23) 


TABLE  6 


r 


nondimensional  hydrodynamic  coefficients 


X* 

u 

xz  — 

0.015020 

Y * = 

V 

- 0.082369 

X*  • 
u 

s=  — 

0.001623 

Y 1 = 

r 

0.026955 

X' 

0 

= _ 

0.000086 

Y'  = 
<P 

0.000086 

x;8 

= _ 

0.002770 

Y f • = 

V 

- 0.035545 

V 

w 

= _ 

0.050138 

Y'.  = 
P 

- 0.000190 

7.' 

0.017455 

y.  = 

0.000400 

4 

z'- 

= __ 

0.031545 

Y'  = 

0.026262 

w 

6R 

z'  • 
q 

= _ 

0.000130 

K'  = 

V 

0.002579 

= - 

0.027695 

K'  = 
P 

- 0.000067 

s 

Mf 

w 

= 

0.009550 

K 1 = 

r 

- 0.000280 

Mo 

= - 

0.001530 

K = 

- 0.001530 

M’ 

q 

= - 

0.011310 

K*  * = 

V 

0.000190 

M*  • 
w 

= - 

0.000146 

K ' • = 
P 

- 0.000092 

M'  • 

q 

= - 

0.001573 

K*-  = 
r 

- 0.000042 

M r 

6 

s 

= - 

0.012797 

K’  = 
6r 

0.000334 

m f 

= 

0.036397 

N’  = 

V 

- 0.018264 

V 

= 

0.000047 

N'  = 
r 

- 0.012497 

X 

i' 

= 

0.001917 

N 1 • = 

V 

- 0.000182 

y 

r 

z 

= 

0.001547 

N 1 • = 
r 

- 0.001531 

It. 

- 

49.33  feet 

N«  = 

- 0.012663 

u 

3 

10.13  ft/sec 

O 

22 

TABLE  7 

NONDIMENSIONAL  TRANSFER  FUNCTION  COEFFICIENTS 


A 

B 

C 

D E 

A, 

Long 

-0.00000902 

-0.00003960 

-0.00003288 

-0.00001028  -0.00000115 

Nf 

s 

0.00000066 

0.00000262 

0.00000128 

0.00000136 

w 1 

M 

s 

0.00000361 

0.00002256 

0.00000996 

0.00000064 

N9 

6 

s 

0.00003291 

0.00004746 

0.00001361 

A 

Lat 

-0.00000003 

-0.00000018 

-0.00000054 

-0.00000183  -0.00000131 

„V 

\ 

-0.00000001 

-0.00000007 

-0.00000015 

-0.00000069 

< 

0.00000013 

0.00000028 

0.00000149 

0.00000233 

r 

N4 

-0.00000013 

-0.00000099 

-0.00000187 

6r 

0/6 

s 

- 1.10~3[. 

7808  s2  + .2311 

s + .0136] 

m 

o 

r-i 
' rH 

1 

0769  s4  + 4 

.5761  s3  + .7799  s2  + .0501  s + .0012] 

Likewise  for  the  v'/6  transfer  function, 

K 


v'/6r  = 


- 1.10~5[.001  s'4  + .007  s'3  + .015  s 


•2  + 


.069  s' 


1.10  5[.003  s'5  + .018  s'4  + .054  s'3  + .183  s'2  + .131  s' 


is  converted  to  dimensional  form  by  the  substitutions  s'  = s X./U  and 

o 

vf  = v/Uq.  Thus,  the  dimensional  transfer  function  is 


v/6r 


.5829  s + .7698  s 


.8460  s5  + 1.003  s4  + 


3 + .3510  s2  + .3383  s 

6272  s3  + .4352  s2  + .0637  s 


The  dimensional  transfer  functions  are  given  in  Table  8 in  the 
polynomial  and  factored  form. 


(Text  Continued  on  Page  25) 


23 


TABLE  8 


DIMENSIONAL  TRANSFER  FUNCTIONS 

1.0. 10~3(0. 7688s3  + 0.6306s2  + 0.0633s  + 0.0014) 

- 1.0. 10~2(0. 5077s4  + 0.4576s3  + 0.0780s2  + 0.0050s  + 0.0001) 

- 0. 1514(s  + 0.0302) (s  + Q.0826)(s  + 0.7074) 

[(s  + . 0592) 2 + (0 . 0219) 2 ] (s  + 0.0811)(s  + 0.7018) 

1.0.10~2(0. 4228s3  + 0.5420s2  + 0.0491s  + 0.000b) 

- 1.0. 10~2(0. 5077s4  + 0.4576s3  + 0.0780s2  + 0.0050s  + 0.0001) 

- Q.8328(s  + 0.0158) (s  + 0.0811)(s  + 1.1849) 

[ (s  + . 0592) 2 + (0.0219)2] (s  + 0.0811)(s  + 0.7018) 

l.o,  10~3(0. 7808s2  + 0.2311s  + 0.0136) 

- 1 .0. 10_2(0. 5077s4  + 0.4576s3  + 0.0780s2  + 0.0050s  + 0.0001) 

_____  - 0. 1538(s  + 0.0811)  (s  + 0.2150) 

[(s  + .0592) 2 + (0.0219)2] (s  + 0.0811) (s  + 0.7018) 

- 1.0.10~4(0. 5829s4  + 0.7698s3  + 0.3510s2  + 0.33835) 

- 1.0.10"4(0. 8459s5  + 1.0033s4  + 0.6272s3  + 0.4352s2  + 0.0637s) 

0 . 6890 (s) [ (s  + . 0516) 2 + (0.6885)2l(s  + 1.2176) 

s[(s  + . 0543) 2 + (0 .6768) 2 ] (s  + 0.1826)(s  + 0.8949) 

-_1 .0.10~4(0. 1460s3  + 0.0658s2  + 0.0725s  + 0.0233) 

- 1.0. 10_4(0. 8459s5  + 1.0033s4  + 0.6272s3  + 0.4352s2  + 0.0637s) 

0 . 1 7 2 5 [ ( s + .0522) 2 + (0.6767)2T(s  + 0.3462) 
s[(s  + . 0543) 2 + (0.6768)2] (s  + 0.1826)(s  + 0.8949) 

- 1.0. 10~4(0. 1461s3  + 0.2340s2  + 0.0910s) 

- 1.0.10_4(0. 8459s5  + 1.0033s4  + 0.6272s3  + 0.4352s2  + 0.0637s) 

0,1727(s)(s  + 0 , 9378) (s  + 0.6638) 

s [ ( s + . 0543) 2 + (0.6768) 2 ] (s  + 0.1826)(s  + 0.8949) 


The  characteristic  equation  for  longitudinal  motions  has  four 
roots.  The  second  order  mode  is  ]ow  frequency  (0.06  rad/sec)  and  very 
highly  damped  (5=0.94).  In  the  forward  velocity  (u)  transfer  functions 
the  zeros  at  0.0826  and  0.7074  very  nearly  cancel  the  poles  at  0.0811 
and  0.7018.  Thus  the  perturbations  in  forward  velocity  will  be  repre- 
sented by  a second  order  response  with  one  overshoot  and  no  undershoot, 
with  the  0.95  value  for  settling  time  equal  to  50.7  secs.  The  response 
in  vertical  velocity  (w)  and  pitch  (9)  are  similar  to  the  u-response, 
the  major  difference  being  the  value  of  the  numerator  zero. 

The  characteristic  equation  for  the  lateral  motions  has  five  roots 
First,  there  is  a free  s in  the  denominator.  Then,  there  is  a second 
order  mode  which  is  of  medium  frequency  (0.68  rad/sec)  and  very  lightly 
damped  (5=0.08).  For  the  side  velocity  (v)  transfer  function,  the  free 
s is  canceled  in  the  numerator  and  the  second  order  mode  nearly  cancels 
the  second  order  mode  in  the  numerator.  Thus,  perturbation  in  side 
velocity  will  have  a first  order  response  mode.  The  characteristic 
mode  will  be  the  root  at  0.1826  with  a 0.95  value  for  settling  time 
equal  to  16.4  seconds. 

For  the  yaw  transfer  function,  the  second  order  mode  again  is 
nearly  canceled  but  the  free  s in  the  denominator  remains.  Thus,  the 
response  will  be  characterized  by  a first  order  integrator.  For  the 
roll  transfer  functions,  the  free  s is  canceled  in  the  numerator  and 
the  second  order  mode  is  the  characteristic  response  mode.  Thus,  the 
motions  in  roll  will  be  oscillatory,  with  a period  of  9.25  seconds  and 
lightly  damped. 

Figures  3 and  4 show  the  time  responses  of  the  vehicle  for  a one- 
degree  step  in  the  control  surface.  These  time  responses  were  computed 
on  a digital  computer  using  the  full  transfer  functions  as  shown  in 
Table  8. 


SUMMARY 


The  linear,  small-perturbation  equations  of  motion  have  been 
derived  by  equating  the  forces  and  moments  acting  on  the  vehicle  to  its 
reactions.  Appropriate  assumptions  are  introduced  to  separate  the 
longitudinal  motions  from  the  lateral  motions.  The  Laplace  transform 
method  of  solution  is  used  to  solve  for  the  independent  variables.  Ex- 
pansion of  the  transfer  functions  polynomial  coefficients  is  completed 
and  the  results  are  presented  in  the  appendices  for  reference.  An 
example  problem  is  provided  to  demonstrate  the  use  of  the  transfer 
functions  and  to  present  typical  values  for  the  characteristic  vehicle 
motions . 
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FIGURE  4.  TIME  RESPONSES  FOR  LATERAL  MOTIONS 


APPENDIX  A 


EXPRESSIONS  FOR  THE  LONGITUDINAL  TRANSFER  FUNCTION  COEFFICIENTS 

The  longitudinal  characteristic  equation  is 

At  = As'4  + Bs'3  + Cs'2  + Ds'  + E 
Long 

where 

A = (m'  - X ' • ) (m'  - Z ' - ) (I ’ - M'*)  - X'-M'-Z'-  - Z'*M'*X'- 
u wy  q wuq  uwq 

- (m'  - Z ’ • ) M ' • X * * - Z’-X'-CI’  - M' • ) - - X'*)Z'*  . 

wuq  u w y q w u q 

B = - (m'  - X'-Xm'  - Z ' • )M'  - Z'(m'  - X'*)(I'  - M'-) 
u w q w u y q 

- X ' (m ' - Z'*)(I'  ~ M'-)  - X ' • M ' • ( Z ' + m')  - M'X'*Z'*  - M'-X'Z'' 

u wy  q wuq  uwq  uwq 

- Z'-M'-X'  - Z’-M'Z'-  - Z'M'-X'-  - (in'  - Z'-)M'-X'  - M'(m'  - Z'-)X'- 

uwq  uwq  uwq  w uq  u w q 

+ M'  *Z  'X'  • + Z'-X'-M'  - Z'X'*(I’  - M ' • ) - Z'-X'd'  - M'*) 
uwq  uwq  uwy  q uw  y q 

- M’-(m'  - X'*)(Z'  -t-  m')  - M' (m’  - X’*)Z’-  + M’-X’Z'*  . 

w uq  w u q wuq 

C = - (m’  - X'  •)  (in'  - Z ' • )M'  + Z'(m'  - X'*)M’  + X'(m'  - Z'*)M' 
u w7  0 w u q u w q 

+ Z'X'(I'  - M'*)  - X'*M'*Zl  - M'X'*(Z'  + m')  - M'*X'(Z'  + in') 

wu  y q wu0  uw  q uw  q 

- Z’*M'-X'  - Z* •M'X'  - Z'M'*X'  - Z'M'X' • - (m'  - Z'*)M'*X! 

u wo  uwq  uwq  uwq  w u9 

- M’(m'  - Z'*)X'  + M'-Z'X'  + Z'M'X' • + Z'*X'*Ml  + Z'X'*M'  + Z'*X'M' 

u w q uwq  wuq  uw0  uwq  uwq 

- X'Z'(I'  - M'*)  - M'*(tn'  - X’*)Z!  - M’(m’  - X’*)(Z’  + m’) 

wuy  q w u 0 w u q 


+ M'*X’(Z'  + ra')  + X'M'Z ' • - M'X' Z ' • . 
wu  q uwq  uwq 


o - w'  - ’t,i> + x;<m'  - z w,m;  - kk%  - - m'lkz'o 

■ wz; + ■’>  - zwj  - z;«'«  - w;  - vm'  - z vx; 

+ + w; + z;x'c,«; + + w;  - - x'szi 


+ M'-X’Zl  + M'X'(Z'  + ra')  . 
w u 0 w u q 


E - - Z'X'M'  - M’X'Z!  - Z'M'X'  + i'M'Xl  + X'Z'Ml  + X'M’Zl 
w u 0 UW0  u v 6 WU0  WU0  u w 0 


The  pitch  response  transfer  function  is 


\ Ae,'2  + Bqs‘  + C9 

0/6S  "A  “ A 

Long  Long 


where 


Ae  = M'  - X--)(m-  - Z'.)  + Z’  X'-M'.  + X’  Z'-M'-  - X'-Z'- 


+ X’  (m'  - Z’^M'.  + V - X'^)  . 

e e 


Be  = - m6  (m’  - x’u)z;  - m6  x;(m’  - z’w>  + za  x’wm;  + z6_x;m': 


u w 6 u 
e 


+ Xl  Z'-M'  + Xl  Z’M’*  - Ml  X';Z'  - M'  X'Z'-  + X'  (id'  - Z'-)M' 

6 u w 6 u w owu  owu  o wu 

e e e e e 

- Xl  Z'M'-  - Zl  M’-X'  + Zl  M’(m'  - X'O  . 

6wu  owu  o w u 

e e e 


c„  = Ml  X'Z'  + Zl  X'M'  + Xl  Z'M'  - MI  X'Z'  - Xl  Z'M’  - Zl  M'X' 
0 6 u w 6 w u 6 u w 6 w u 6 wu  5wu 

e e e 


The  vertical  velocity  transfer  function  is 


A s'  + B s'  + C s'  + D 
w w w w 


6 A. 
s Long 


where 


A = Zi  (m'  - X’*)(I’  - M'*)  + Xi  M' • Z ' • + Mi  Z'*X'*  - Zi  M'*X’* 
w 6 u y q 6uq  ftuq  ftuq 

e yne  ene 

+ X!  Z * - (I  * - M ' • ) + Mi  (m'  - X'*)Z'*  . 

ft  u y q ft  u q 

e e n 

B = - Zi  - X ' • )M ' - Zi  X'(I'  - M’*)  + Xi  M’Z'*  + Xi  M'-(Z'  + m') 
w 6 u q ft  u y q 6uq  ftuq 

e e e e 

+ Mi  Z'-X'  + Mi  Z'X'-  - Zi  M'-X'  - Zi  M'X'  • - Xi  Z'-M' 
ft  uq  ftuq  ft  uq  ftuq  6 uq 
e e e en  e 

+ Xi  Z * (I ' - M'*)  + Mi  (m'  - X’*)(Z'  + m')  - Mi  X'Z'*  . 
ft  u y q 6 u q ftuq 

e e M e 

C = - Zi  (m'  - X' • )M'  + Zi  X’M'  + Xi  M'*Zi  + Xi  M'(Z'  + m')  + M'  Z'*X' 

w ft  u 6 ftuq  6u6  ftuq  Su0 

e e e e e 

+ Mi  Z'Z’  - Z'  M'*Xi  - Zi  M'X'  - Xi  Z'*Mi  - Xi  Z'M'  + Mi  (m'  - X'*)Zi 

6uq  5 u 0 6 u q 6u0  6uq  6 u0 

e e e e e e 

- M’  X’(Z’  + ra* ) . 

6 a q 
e 

D = Z\  X’M’  + X*  M’Z’  + M’  Z’X’  - Zi  M’X’  - X’  Z’M’  - M’  X’Z’  . 

w 5u°  6 u 9 5u0  5u0  ftuo  6u0 

e e e e e e 

The  forward  speed  transfer  function  is 


A s'3  + B s'2  + C s'  + D 


u u 


6 A 

s Long 


whe  re 


A = Xi  (m’  - Z'-)(I'  - M'-)  + Mi  X’*Z'-  + Zi  M'-X’*  + M'  (in'  - Z'OX'* 
u 5 w y q ftwq  ftwq  wq 

e e e e 


+ Zi  X*  * (I  * - M'-)  - Xi  M'-Z'*  . 
ft  w y q ftwq 
e e 


mmm 


B - - Xl  (m* 
u 6 

- Z'-)M'  - XI  ( 1 ' - 
w q 6 y 

M1  • )Z  * + 
q w 

Mi  X' • (Z f + raf ) + Ml  X*  Zf • 
6 w q 6 w q 

e 

e 

e 

e 

+ Z!  M 1 • X 1 

+ Zl  M ’ X ' • + Ml  (m' 

- Z'-)Xf 

- M!  ZfX'*  - Z 

: x'-m* 

6 w q 

e 

6 w q 6 

e e 

w q 

6 w q 

e 

6 w q 

e 

+ Zl  X'(I'  - M'*)  - XI  M' • (Z ' + m')  - XI  M'Z'*  . 
6 w y q 6 w q 6 w q 


APPENDIX  B 


EXPRESSIONS  FOR  THE  LATERAL  TRANSFER  FUNCTION  COEFFICIENTS 

The  lateral  characteristic  equation  is 

ALa(;  = s'(As'4  + Bs'3  + Cs'2  + Ds'  + E) 


where 


A = (m-  - Y'.)(i;  - N'.)(r  - K-.)  + N'.(-  1^  - K’^Y’. 

+ Y’rK'v(-  Kz  ~ N'p>  ~ N'vY’r(K  ~ K’p>  “ K'v(i;  ' N'^Y'p 


- (m’  -YV<-  -»V("  • 


B = - (m’  - Y 1 • ) (I ’ - N’-)K’  - N’(m'  - Y'-)(I'  - K'-) 
v z r'  p r v'v  x p' 

- Y'(i;  - N'-)(I'  - K’;)  + N ’ • (-  I'  - K'.)Y'  - N'-R’Y’- 

vz  rx  p v xz  r p vrp 

+ N’(-  I'  - K'*)Y'«  - Y ' »K  ' *N  ' + Y ' • K'  (-  I'  - N'O 
v xz  r p r v p r v xz  p’ 

- (m-  _ y;)K’.(-  rz  - N’.)  + N'-Y’.K’  + N'.(m'  - Yp(i;  - K'.) 

- Vr(i;  - K’p>  - K'v(Iz  “ N'r)Yp  + K';N;Y’p  - " N’ W’ 

+ (m'  - Y ' • ) (-  I'  - K’;)N'  + (m'  - Y'-)K'(-  I'  - N'*) 
v xz  r p v r xz  p 

+ Y ' (-  I*  - K' • ) (-  I’  - N'*)  . 
v xz  r xz  p 


B-l 


C =•  N'(m'  - Y ' • )K'  - (m'  - Y'-)(I'  " N'*)K!  + Y * (I ' - N'*)K' 
r v p v z r <p  v z rp 

+ N ' Y ' ( I ' - K'*)  + N ' • (-  I'  - K'*)Y!  - N ' ‘ K ' Y ' + N'(-  I'  - K'-)Y' 
r v x p v xz  r v r p v xz  r p 


- N'K'Y’ 

'•  - Y’-K'-N! 

v r 

p r v 

- (m'  - 

Y')K'(-  I' 

r v xz 

+ N' (m' 

V 

- Y ' ) ( I ' - K 
r x 

r v p 


r v p 


r p v r p 


v r p z r v p 


+ K'N'Y'-  + (m'  - Y ' • ) (-  I’  - K*  * )N ' - (m’  - Y’-)K'N' 
vrp  v xz  r<p  vrp 


- Y'(-  1’  - K' • )N ' - K'Y'(-  I'  - N'-)  . 

v xz  r p r v xz  p 


D = N'(ra'  - Y'-)K’  + Y ' (I ' - N'-)K!  - N'Y'K'  - N'-K'Y! 
r v <p  v z r <J>  rvp  v r <f> 

+ N'(_  x»  _ R ' • ) Y ' - N'K'Y'  - Y'-K'N'  + (m*  - Y')K'-N' 
v xz  r <t>  vrp  r v 4>  r v 

+ (m*  - Y')K'N’  - N ' • (m'  - Y')K!  + N'Y'-K!  - N'(ra'  - Y')K'  + K';N'Y' 
r vp  v r<t>  v r <p  v rp  vr<j 

- ( 1 ' - N'-)K'Y!  + K'N’Y'  - (m’  - Y'-)K'N!  - Y'(-  I'  - K';)N! 
z r v <p  vrp  v r <p  v xz  r <?> 


+ K'Y'N'  . 
rvp 


E’  = - N'Y'K' 
r v <p 


N'K’Y'  + (in'  - Y ’ ) K 'N  ' - N'  (m' 
v r <t>  r v v 


Y')K'  + K'N'Y’  + K'Y'N' 
r d>  vr<f>  rvd 


'fhe  sideslip  transfer  function  is 


s'<v3  + + v 

5R  ^Lat  ALat 


where 


A,  = Y!  (I'  - N ' • ) (I ' - K' • ) - Nl  (-  I'  - K'’)Y'*  - K!  Y'*(-  I'  - N'*) 
B 6 z rx  p 6 xz  r p 6 r xz  p 

r v r K r r 

+ N?  Y ' • (I ' - K'«)  + Kl  (I’  - N'-)Y'>  - Yl  (-  I'  - K' • ) (-  I'  - N'*) 
6 r x p 6 z r p 6 xz  r xz  p 
r r r r r 

B„  = - YI  (I'  - N ' • )K'  - Y!  N'(I'  ~ K ' • ) - Nl  (-  I'  - K'-)Y'  + Nl  K'Y' • 
3 6 z r p 6 r x p 6 xz  r p 6 r p 

r r r r r 

+ K'.  Y ' *N ' + KI  (m’  - Y ' ) (-  I’  - N’-)  - N!  Y'*K' 

6rp  6 r xz  p orp 

r r r r r r 

- (m'  - Y')NI  (I'  - K'*)  + K!  (I'  - N'OY'  - K!  N ' Y 1 • 
r 6 x p o z r p 6 r p 
r K r r 

+ Y!  (-  I*  - K’-)N'  + Y ' K'(-  I'  - N'-)  . 

5 xz  r p 6 rx  xz  p 

r r 

C0  = - Y!  (I'  - N ' • )K'  + Y!  N'K'  - N'  (-  I'  - K'*)Y'  + N'  K'Y' 

6 5 z r <p  6rp  6 xz  r 5rp 

r y r r r 

+ K!  Y'*N!  - KI  (m'  - Y')N'  - N!  Y'-K!  + (m'  - Y')N!  K' 

6ri>  6 r p o r <j>  rop 

r r r r 

+ K’  (I*  - N’*)Y!  - N'K’  Y'  + YI  (-  I'  - K'-)N!  - Y'  K'N ' . 

6 z r d>  r6p  6 xz  r ♦ 6rp 

r r r r 

Dn  = Y'  N'K!  + Nl  K'Y!  - Ki  (m'  - Y')N'  + (m'  - Y')N’  K!  - N'K'  Y! 
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